CELLmicrocosmos 2.2: advancements and applications in modeling of three-dimensional PDB membranes

Sommer B, Dingersen T, Schneider S, Rubert S, Gamroth C. CELLmicrocosmos 2.2: advancements and applications in modeling of three-dimensional PDB membranes (Conference Abstract). In: Journal of Cheminformatics. Journal of Cheminformatics. Vol 2(Suppl 1):O21. Springer Science and Business Media LLC; 2010

Block of NMDA receptor channels by endogenous neurosteroids: implications for the agonist induced conformational states of the channel vestibule

N-methyl-D-aspartate receptors (NMDARs) mediate synaptic plasticity, and their dysfunction is implicated in multiple brain disorders. NMDARs can be allosterically modulated by numerous compounds, including endogenous neurosteroid pregnanolone sulfate. Here, we identify the molecular basis of the use-dependent and voltage-independent inhibitory effect of neurosteroids on NMDAR responses. The site of action is located at the extracellular vestibule of the receptor's ion channel pore and is accessible after receptor activation. Mutations in the extracellular vestibule in the SYTANLAAF motif disrupt the inhibitory effect of negatively charged steroids. In contrast, positively charged steroids inhibit mutated NMDAR responses in a voltage-dependent manner. These results, in combination with molecular modeling, characterize structure details of the open configuration of the NMDAR channel. Our results provide a unique opportunity for the development of new therapeutic neurosteroid-based ligands to treat diseases associated with dysfunction of the glutamate system

Use of molecular modelling to probe the mechanism of the nucleoside transporter NupG.

Nucleosides play key roles in biology as precursors for salvage pathways of nucleotide synthesis. Prokaryotes import nucleosides across the cytoplasmic membrane by proton- or sodium-driven transporters belonging to the Concentrative Nucleoside Transporter (CNT) family or the Nucleoside:H(+) Symporter (NHS) family of the Major Facilitator Superfamily. The high resolution structure of a CNT from Vibrio cholerae has recently been determined, but no similar structural information is available for the NHS family. To gain a better understanding of the molecular mechanism of nucleoside transport, in the present study the structures of two conformations of the archetypical NHS transporter NupG from Escherichia coli were modelled on the inward- and outward-facing conformations of the lactose transporter LacY from E. coli, a member of the Oligosaccharide:H(+) Symporter (OHS) family. Sequence alignment of these distantly related proteins (∼ 10% sequence identity), was facilitated by comparison of the patterns of residue conservation within the NHS and OHS families. Despite the low sequence similarity, the accessibilities of endogenous and introduced cysteine residues to thiol reagents were found to be consistent with the predictions of the models, supporting their validity. For example C358, located within the predicted nucleoside binding site, was shown to be responsible for the sensitivity of NupG to inhibition by p-chloromercuribenzene sulphonate. Functional analysis of mutants in residues predicted by the models to be involved in the translocation mechanism, including Q261, E264 and N228, supported the hypothesis that they play important roles, and suggested that the transport mechanisms of NupG and LacY, while different, share common features

Pharmacophore elements of the TIPP class of delta opioid receptor antagonists

A series of tri-and tetrapeptides sharing the amino-terminal dipeptide unit Tyr-Tic, found in the high-affinity delta opioid receptor antagonist Tyr-Tic-Phe-Phe (TIPP), was prepared and evaluated in receptor binding assays to explore the role(s) of the phenylalanine residues in positions 3 and 4. It was found that aromaticity of residues 3 and 4 is not required for high affinity, a lipophilic side chain in either location being sufficient, as evidenced by the high delta receptor binding affinities observed for the tetrapeptide Tyr-Tic-Ala-Leu and the tripeptide Tyr-Tic-Leu. These results support the suggestion of Temussi et al. [Biochem. Biophys. Res. Commun., 198 (1994) 933] that the aromatic side chain of the Tic residue corresponds to the aromatic side chain found in residues 3 or 4 in other delta-selective peptide series.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/43172/1/10989_2004_Article_BF00126275.pd

Structure-based statistical analysis of transmembrane helices

Recent advances in determination of the high-resolution structure of membrane proteins now enable analysis of the main features of amino acids in transmembrane (TM) segments in comparison with amino acids in water-soluble helices. In this work, we conducted a large-scale analysis of the prevalent locations of amino acids by using a data set of 170 structures of integral membrane proteins obtained from the MPtopo database and 930 structures of water-soluble helical proteins obtained from the protein data bank. Large hydrophobic amino acids (Leu, Val, Ile, and Phe) plus Gly were clearly prevalent in TM helices whereas polar amino acids (Glu, Lys, Asp, Arg, and Gln) were less frequent in this type of helix. The distribution of amino acids along TM helices was also examined. As expected, hydrophobic and slightly polar amino acids are commonly found in the hydrophobic core of the membrane whereas aromatic (Trp and Tyr), Pro, and the hydrophilic amino acids (Asn, His, and Gln) occur more frequently in the interface regions. Charged amino acids are also statistically prevalent outside the hydrophobic core of the membrane, and whereas acidic amino acids are frequently found at both cytoplasmic and extra-cytoplasmic interfaces, basic amino acids cluster at the cytoplasmic interface. These results strongly support the experimentally demonstrated biased distribution of positively charged amino acids (that is, the so-called the positive-inside rule) with structural data

Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics Simulations

Molecular dynamics simulations of membrane proteins have provided deeper insights into their functions and interactions with surrounding environments at the atomic level. However, compared to solvation of globular proteins, building a realistic protein/membrane complex is still challenging and requires considerable experience with simulation software. Membrane Builder in the CHARMM-GUI website (http://www.charmm-gui.org) helps users to build such a complex system using a web browser with a graphical user interface. Through a generalized and automated building process including system size determination as well as generation of lipid bilayer, pore water, bulk water, and ions, a realistic membrane system with virtually any kinds and shapes of membrane proteins can be generated in 5 minutes to 2 hours depending on the system size. Default values that were elaborated and tested extensively are given in each step to provide reasonable options and starting points for both non-expert and expert users. The efficacy of Membrane Builder is illustrated by its applications to 12 transmembrane and 3 interfacial membrane proteins, whose fully equilibrated systems with three different types of lipid molecules (DMPC, DPPC, and POPC) and two types of system shapes (rectangular and hexagonal) are freely available on the CHARMM-GUI website. One of the most significant advantages of using the web environment is that, if a problem is found, users can go back and re-generate the whole system again before quitting the browser. Therefore, Membrane Builder provides the intuitive and easy way to build and simulate the biologically important membrane system

Substrate binding and translocation of the serotonin transporter studied by docking and molecular dynamics simulations

The serotonin (5-HT) transporter (SERT) plays an important role in the termination of 5-HT-mediated neurotransmission by transporting 5-HT away from the synaptic cleft and into the presynaptic neuron. In addition, SERT is the main target for antidepressant drugs, including the selective serotonin reuptake inhibitors (SSRIs). The three-dimensional (3D) structure of SERT has not yet been determined, and little is known about the molecular mechanisms of substrate binding and transport, though such information is very important for the development of new antidepressant drugs. In this study, a homology model of SERT was constructed based on the 3D structure of a prokaryotic homologous leucine transporter (LeuT) (PDB id: 2A65). Eleven tryptamine derivates (including 5-HT) and the SSRI (S)-citalopram were docked into the putative substrate binding site, and two possible binding modes of the ligands were found. To study the conformational effect that ligand binding may have on SERT, two SERT–5-HT and two SERT–(S)-citalopram complexes, as well as the SERT apo structure, were embedded in POPC lipid bilayers and comparative molecular dynamics (MD) simulations were performed. Our results show that 5-HT in the SERT–5-HTB complex induced larger conformational changes in the cytoplasmic parts of the transmembrane helices of SERT than any of the other ligands. Based on these results, we suggest that the formation and breakage of ionic interactions with amino acids in transmembrane helices 6 and 8 and intracellular loop 1 may be of importance for substrate translocation

Structure and dynamics of the active Gs-coupled human secretin receptor

The class B secretin GPCR (SecR) has broad physiological effects, with target potential for treatment of metabolic and cardiovascular disease. Molecular understanding of SecR binding and activation is important for its therapeutic exploitation. We combined cryo-electron microscopy, molecular dynamics, and biochemical cross-linking to determine a 2.3 Å structure, and interrogate dynamics, of secretin bound to the SecR:Gs complex. SecR exhibited a unique organization of its extracellular domain (ECD) relative to its 7-transmembrane (TM) core, forming more extended interactions than other family members. Numerous polar interactions formed between secretin and the receptor extracellular loops (ECLs) and TM helices. Cysteine-cross-linking, cryo-electron microscopy multivariate analysis and molecular dynamics simulations revealed that interactions between peptide and receptor were dynamic, and suggested a model for initial peptide engagement where early interactions between the far N-terminus of the peptide and SecR ECL2 likely occur following initial binding of the peptide C-terminus to the ECD

A membrane-inserted structural model of the yeast mitofusin Fzo1

Mitofusins are large transmembrane GTPases of the dynamin-related protein family, and are required for the tethering and fusion of mitochondrial outer membranes. Their full-length structures remain unknown, which is a limiting factor in the study of outer membrane fusion. We investigated the structure and dynamics of the yeast mitofusin Fzo1 through a hybrid computational and experimental approach, combining molecular modelling and all-atom molecular dynamics simulations in a lipid bilayer with site-directed mutagenesis and in vivo functional assays. The predicted architecture of Fzo1 improves upon the current domain annotation, with a precise description of the helical spans linked by flexible hinges, which are likely of functional significance. In vivo site-directed mutagenesis validates salient aspects of this model, notably, the long-distance contacts and residues participating in hinges. GDP is predicted to interact with Fzo1 through the G1 and G4 motifs of the GTPase domain. The model reveals structural determinants critical for protein function, including regions that may be involved in GTPase domain-dependent rearrangements

Transmembrane protein topology prediction using support vector machines

Background: Alpha-helical transmembrane (TM) proteins are involved in a wide range of important biological processes such as cell signaling, transport of membrane-impermeable molecules, cell-cell communication, cell recognition and cell adhesion. Many are also prime drug targets, and it has been estimated that more than half of all drugs currently on the market target membrane proteins. However, due to the experimental difficulties involved in obtaining high quality crystals, this class of protein is severely under-represented in structural databases. In the absence of structural data, sequence-based prediction methods allow TM protein topology to be investigated.Results: We present a support vector machine-based (SVM) TM protein topology predictor that integrates both signal peptide and re-entrant helix prediction, benchmarked with full cross-validation on a novel data set of 131 sequences with known crystal structures. The method achieves topology prediction accuracy of 89%, while signal peptides and re-entrant helices are predicted with 93% and 44% accuracy respectively. An additional SVM trained to discriminate between globular and TM proteins detected zero false positives, with a low false negative rate of 0.4%. We present the results of applying these tools to a number of complete genomes. Source code, data sets and a web server are freely available from http://bioinf.cs.ucl.ac.uk/psipred/.Conclusion: The high accuracy of TM topology prediction which includes detection of both signal peptides and re-entrant helices, combined with the ability to effectively discriminate between TM and globular proteins, make this method ideally suited to whole genome annotation of alpha-helical transmembrane proteins